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We consider a simple quantum system subjected to a classical random force. Under certain 
conditions it is shown that the noise-averaged Wigner function of the system follows an integro- 
differential stochastic Liouville equation. In the simple case of polynomial noise-couplings this 
equation reduces to a generalized Fokker-Planck form. With nonlinear noise injection new "quantum 
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Abstract 



' diffusion" terms arise that have no counterpart in the classical case. Two special examples that are 

I not of a Fokker-Planck form are discussed: the first with a localized noise source and the other with 

, a spatially modulated noise source. 
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^1 1. Stochastic Liouville Equations 

Stochastic equations have long been used in physics to model various phenomena. 
^ ! Brownian motion, spin relaxation, and critical dynamics may be cited as obvious 
examples. At a formal level there are two ways to set up such equations (1) as exact 
equations or (2) as part of a phenomenological description [Q. In either case one 
typically encounters equations that are nonlocal in time and involve stochastic forcing 
terms usually called "noise." Such Langevin equations exist at both the classical and 
quantum levels. As expected the situation is more complicated in the latter case; 
while in classical problems it is often possible to approximate the "noise" as being 
Gaussian and white, and further to replace a nonlocal kernel by one local in time 
(the Markov approximation), such simplifications do not easily obtain in quantum 
mechanics. Nevertheless, simple approximate approaches are valuable in that they 
often capture some essential physics, or even make some technical point, with less 
calculational clutter when compared to a more comprehensive or refined method of 
attack. The work outlined here is in this spirit. It owes much to Kubo's study of 



the stochastic Liouville equation |^ and a presentation of it given by Zwanzig 0. 
Different aspects of this work have been considered in detail elsewhere ^j. Nonlinear 
couplings to an oscillator environment have been studied in the independent oscillator 
model in Ref. where quantum diffusion has also been shown to exist. 

In this paper, all quantum calculations will be done in the Wigner framework of 
quantum mechanics. Partly this is because quantum distribution functions defined 
on a mock phase space can be easily compared to their classical counterparts. Fur- 
thermore, in the models that will be discussed, stochastic Liouville equations written 
in terms of the Wigner function will be obtained directly from the stochastic Hamil- 
tonian. This enables us to bypass the somewhat delicate question of how to derive 
quantum Fokker-Planck equations starting from Langevin equations for quantum op- 
erators. A nice feature of the phase space approach is that the quantum derivation 
of the stochastic Liouville equation closely parallels the classical derivation; there is 
no need to invoke path integrals. Finally, this approach also enables us to discuss the 
singular nature of the ^ — > limit for both the systematic and the diffusive terms in 
the stochastic Liouville equation. 

We begin with the Hamiltonian (a generalization of the randomly forced oscillator 
considered earlier by Merzbacher [Q): 

H = ^ + V{x)~F{t)g{x), (1) 

where p, x are the dynamical variables characterizing the motion of the system. The 
functions V{x) and g{x) are assumed to be differentiable. F{t) is an external pertur- 
bation that is taken to be Gaussian, white noise, i.e., {F(t))j^ = 0, and 

{F{t,)F{t2))^ = 2B{t,)S{t,-h), (2) 

with the usual restrictions on the higher moments. The ( )^ denotes an average 
over the realizations of F. The delta function in is supposed never to be exactly 
realized, but is treated just as an idealization of a sharply peaked, symmetric function. 
This corresponds to interpreting the noise in the sense of Stratonovich 0. 

One way to write the equations of motion is to use the Liouville equation for the 
phase space distribution function. We introduce the distribution function fci{x,p]t), 
which satisfies the probability flux conservation equation (Liouville's theorem), 
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the right hand side of @ defining the Liouville operator Lci- 

Following Kubo's analysis 0] applied to the Hamiltonian (|I]), we proceed to derive 
the noise-averaged stochastic Liouville equation. With Lq the Liouville operator 
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corresponding to the systematic part of the evolution, we obtain 
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a Fokker-Planck equation for the noise-averaged distribution function. Since fci is a 
phase space distribution function, (Q) is a two-variable Fokker-Planck, or Kramers, 
equation. In the absence of noise it reduces to the usual Liouville equation. We 
observe that whatever V{x) and g{x) may be, ifciit))]\j will always satisfy a Fokker- 
Planck equation. This will not be true in the quantum case, to which we now proceed. 

As in the classical case we will work with the stochastic Hamiltonian ([l|). Be- 
cause of the noise, this Hamiltonian will evolve pure states to mixed states. Thus 
it is appropriate to study not the time dependent Schrodinger equation but rather 
the quantum Liouville equation for the density matrix, given here in the coordinate 
representation, 

ih—p{xi, X2) = {H{xi) - H{x2)*] p{xi,X2). (5) 

We wish to write (^) in the Wigner formalism of quantum mechanics [|l^ and then 
to noise average just as in the classical case. This derivation is given in the first and 
third papers of Ref. and here we quote only the final result: 

{fwiX,k;t))j^ = -LsysifwiX,k;t))j^- / dp {fwiX, k + p;t)) KsiX,p;t), 

J —00 

(6) 
where 
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and Lsys is the systematic quantum Liouville operator. When g{X) can be profitably 
Taylor expanded, the above equation can be written as 
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2. Quantum Diffusion 

The conditions under which @ will reduce to a Fokker-Planck form are when 
both V{X) and giX) are of the form Ax + Bx^. In this case the quantum Liouville 
equation reduces to the classical one. The difference between the two then lies not in 
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the dynamical equation, but in the different constraints imposed on the initial value 

of the respective distribution functions. 

We now study different choices for g{X). If g{X) = AX, with A a constant, then 
= A'^d^/dk'^, a conventional diffusion term. This gives rise to the simple model 



equation often employed in studies of quantum decoherence |]1T 
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Consider now the case, g{X) = AX + \eX^, where 



V = {A + eX2)2^ - -{A + eX')eh^^ + —e^h"^. (10) 
^ ' dk^ 6 ^ ' dk^ 144 dk^ ^ ' 

Notice the appearance of the purely quantum mechanical, higher even derivative 
"diffusion" terms. The classical limit ^ is singular not only for the systematic 



quantum Liouville operator Lsys ITHITH but also for the stochastic terms arising 



from quantum diffusion. It is easy to see that all the quantum diffusive terms, when 
acting on "fast" (cf. Refs. |T^|T^) pieces ~ exp{ikX/h) of a Wigner function, are of 



0{l/?i^). The highest order quantum diffusion term dominates at large distances and 
always acts to increase the linear entropy 1 — / dXdkp [|| . The effect of the quantum 
diffusion terms with regard to decoherence is to reduce the decoherence time at large 
length scales [|| . 
3. Two Illustrative Examples 

As we have seen, the stochastic quantum Liouville equation written in terms of the 
Wigner distribution function is in general a complicated integro-differential equation. 
If the coupling to the noise is through a polynomial in the system variable, then 
this equation can truncate to a finite order partial differential equation. However, 
there are cases of physical interest where the coupling to the noise cannot be reduced 
to such a form. We will now exhibit two such cases, coupling the system (1) to a 
localized noise source, and (2) to a spatially modulated noise source. The first case is 
of interest in quantum tunneling through a stochastic barrier while the second applies 
to the noise in a microwave cavity. More details can be found in the third paper of 
Ref. §. 

A localized noise source can be modeled by setting g{X) = A exp(— eX^/2). In 
this case, 

— {fw{X,k;t))^ = -Lsys{fw{X,k;t))j^ ^Ts / dp {fw{X, k;t)) ^ 
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A spatially modulated noise source, g{X) = asm{i3X/h), leads to 
d 

{fw{X,k;t))j^ = —Lsys {fw{X,k;t))j^ 
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4 {fw{X, k - (3- 1))^ - I {fwiX, k + (5- t))j. 



(12) 



Eqs. ([Tl| ) and (|T^) are not of a classical form: the corresponding classical equations 
result from keeping only the first term of a derivative expansion of fw in these equa- 
tions (the quantum equations may be viewed as resulting from a resummation of all 
terms in such an expansion). 
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